arXiv:1506.07470v4 [math.OA] 12 Jan 2017 


DYNAMICAL DEFORMATION OF TOROIDAL MATRIX 

VARIETIES 

FREDY VIDES 


Abstract. In this document we study the local connectivity of the sets whose 
elements are m-tuples of pairwise commuting normal matrix contractions. 
Given e > 0, we prove that there is <5 > 0 such that for any two m-tuples 
of pairwise commuting normal matrix contractions X := {Xi,.. ., Xm) and 
X := (Xi,..., Xm) that are 5-close with respect to some suitable distance 9 
in ggj g m-tuple of matrix paths (homotopies) connecting 

X to X relative to the intersection of some e, 9-neighborhood of X with the 
set of m-tuples of pairwise commuting normal matrix contractions. One of the 
key features of these matrix homotopies is that 5 can be chosen independent 
of n. 

Some connections with topology and numerical matrix analysis will be out¬ 
lined as well. 


1. Introduction 

In this document we study local deformation properties of what we call toroidal 
matrix varieties. Given a m-tuple of u x n matrices X S (^£nxn^m ^ suitable 
distance S in induced by the spectral norm (that will be defined below), 

let us write Ag(X, r) to denote the set {Y s (C”^")™|6(X, Y) < r}. We call 
Ag(X, r) a r, 9-neighborhood of X. 

Given two matrices X,Y G let us write X Y to indicate that there is 

a homotopy between X and Y, in other words, there is function 7 S C([0,1], C"^") 
such that 7(0) = X and 7(1) = Y (where C{X, Y) denotes the algebra of continuous 
functions on X that take values in Y). Given e, 5 > 0, any two m-tuples of pairwise 
commuting normal matrix contractions X := {Xi, ..., Xm) and X := {Xi, ..., Xm) 
such that 9(X, X) < e, let us write X X to denote the homotopy induced by the 
homotopies of their components. We study the existence of homotopies X X 
relative to the intersection of some e, 9-neighborhood of X with the set of m-tuples 
of pairwise commuting normal matrix contractions. 

We can think of each homotopy Xj Xj relative to Ag(Aj, e), as a noncommu- 
tative analogy of the family of curves (links) connecting the point sets determined 
by the spectra cr{Xj) and <j{Xj) relative to the embedding of x in [—1,1]^ x 
(where the embedding is induced by the mapping —>■ i?(0,1) C [—1,1]^). In a 
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similar way we can interpret the induced homotopies X X as noncommuta- 
tive analogies of the links connecting the joint spectra cr(X) and o'(X) relative to 
X Tk 

In this study we will use the (7*-algebraic technology developed in [iniiiiiiiiiis] 
to derive a connectivity technique that can be used to study and solve problems 
that appear in clustered matrix approximation (in the sense of |28p and matrix 
dilation/compression problems (in the sense of [121 [T4] and [TJ §9 and §10]). 

In [22] we found some connections between the map 'I' := Ad[bk] described by 
LEU in ]J2] and the local homeomorphims that play a key role in the Kirby Torus 
Trick (introduced by R. Kirby in |19j l by moving sets of points in homeomorphic 
copies T™ — T’" of T™ just a little bit, this analogy together with L I3.il provides 
a connection with topologically controlled linear algebra (in the sense of m), 
which can be roughly described as the study of the relations between matrix sets 
and smooth manifolds that is performed by implementing techniques from geomet¬ 
ric topology in the study of matrix approximation problems. The search for the 
previously mentioned analogies and connections was motivated by a question raised 
by M. H. Freedman regarding to the role played by the Kirby Torus Trick in linear 
algebra and matrix approximation theory. 

In this document we build on the techniques developed in [22] and [29] to study 
the analytic local connectivity properties of matrix representations of certain uni¬ 
versal C*-algebras. The deformation technology constructed for this purpose can be 
adapted using the techniques presented in ^to study the local behavior of approx¬ 
imate/numerical solutions to matrix equations in the sense of [HJIHIISI |20j [HI [18], 
that can be interpreted as local matrix homotopies that play the role of continuous 
analogies (in the sense of m of the discrete/numerical procedures corresponding 
to spectral refinement algorithms (in the sense of mm that can be used to re¬ 
fine the outputs of approximate simultaneous block diagonalization/transformation 
algorithms (in the sense of [23] [24] ). 

The applications of the previously described refinement methods include singular 
saddle point problems that arise in computational fluid dynamics, mixed finite ele¬ 
ment approximation of elliptic partial differential equations, optimization, optimal 
control, weighted least-squares problems, electronic networks, computer graphics 
and others in scientific computing and engineering (in the sense of [20] 1. together 
with mathematical modelling of photonic crystals (in the sense of m) in physics. 
The results presented in ^can also be used to compute and to study the approxi¬ 
mate and exact simultaneous block diagonalization of matrix m-tuples in the sense 
of [231121], problems of this type appear in biomathematics, image processing and 
applied spectral graph theory. 

Building on some of the ideas developed by M. A. Rieffel in [261127], by H. Lin in 
[21] and by P. Friis and M. Rordam in [10], we proved L I3.il L I3.2l and L I3.31 these 
results together with L I2.il provide us with the matrix approximation technology 
that we use to prove the main results. The main results T I3.il C I3.1l and T I3.2l are 
presented in ^ and some future directions are outlined in El 

2. Preliminaries and Notation 

Let {X, d) be a metric space. Then we say that C A is a d-dense if for all 
X G X there exists x G X such that d{x, x) < 6. Given two compact subsets X, Y of 
the complex plane, we will write d_f/(A, F) to denote the Hausdorff distance between 
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X and Y that is defined by the expression dniXjY) := max{supa;eA: infj/gy |a; — 
y\, supy£Y \y—x\}. Let us denote by S the function S : —>■ 

IRq,(S,T) i-a maxj \\Sj — Tj\\. From here on || • || denotes the operator/spectral 
norm. Given two topological objects, we use the expression X ~ F to denote 
a homeomorphism between them (a continuous function between the topological 
objects that has a continuous inverse function). We will write D"* and T™ to denote 
the m-dimensional closed unit disk and the m-dimensional torus respectively. We 
will write B{xo, r) to denote the closed ball {a; € C | |a; — xq] < r}. 

We will write Mn to denote the set C"^” of n x n complex matrices, the sym¬ 
bols In and On will be used to denote the identity matrix and the zero matrix in 
Mn respectively. Given a matrix A G Mn, we write A* to denote the conjugate 
transpose A^ of A. A matrix X G Mn is said to be normal if XX* = X*X, a 
matrix H G Mn is said to be hermitian if H* = H and a matrix U G Mn such that 
U*U = UU* = In is called unitary. 

A hermitian matrix P G Mn such that P = is called a projector. Given two 
projectors P and Q, if PQ = QP — 0„ we say that P and Q are orthogonal. By an 
orthogonal partition of unity in Mn, we mean a finite set of orthogonal projectors 
{Pj} in Mn such that Pj = 1„. We will omit the explicit reference to Mn when 
it is clear from the context. In this document by a matrix contraction we mean a 
matrix X in Mn for some n G Z+ such that ||A|| < 1. 

Given any two matrices X,Yg Mn we will write [A, Y] and Ad[A](F) to denote 
the operations [A, A] := AA - AA and Ad[A](A) := AAA*. 

Given a matrix A G Mn, we write a{X) to denote the set {A G C| det(A —Al„) = 
0} of eigenvalues of A, the set cr(A) is called the spectrum of A. Given a compact set 
X c C and a subset S C M„, let us dote by S'(X) the set 5'(X) := {A G S'|a-(S') C X}, 
for instance, the expression A/'(n)(D^) is used to denote the set normal contractions 
in Mn- 

Definition 2.1. Given e > 0 and a matrix X G Mn, we write cre(A) to denote the 
e-Pseudospectrum of X which is the set defined by the following relations. 

(Te(A) := {A G C I A G cr(A-I-A), for some if G M„ with ||if|| < e} 

= {A G C I II Av — Av|| < e, for some v G C" with ||v|| = 1} 

Definition 2.2 (Semialgebraic Matrix Varieties). Given J G Z+, a system of J 
polynomials pi,... ,pj G LI^jv) = C (a;i,... ,xjv} in N NG-variables xi,... ,xn G 
n^jv) dnd a real number s > 0, a particular matrix representation of the noncom- 
mutative semialgebraic set n{pi, ■ ■ ■ ,Pj) described by 
(2.1) 

Ze,n{pi,---,pj) := { Ai,..., Aa, G M„ I \\pj{Xi,...,XN)\\ < e,l < j < J} , 

will be called a e, n-semialgebraic matrix variety (e,n-SMV), if e = 0 we can 
refer to the set as a matrix variety. 

Example 2.1. Given any integer n > 1, let us set N := diag [n, n — 1,..., 1], we 
will have that the set An := {A G M„|[N,A] = 0} zs o matrix variety. If for 
some S > 0, we set now An ,5 := {A G ilin| || [N, A] || < <5}, the set An ,5 is a matrix 
semialgebraic variety. 

Example 2.2. The subset of Mn described by, 

(2.2 ) A3:=|(Ai,A2,A3)gM3 


^lAi — X2B1 — Cl, 1 
GI2A3 — X2B2 = C2 j ’ 
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where Ai,A 2 ,Bi,B 2 ^Ci and C 3 are some fixed but arbitrary matrix contractions 
in Mn, is an algebraic matrix variety. 

Example 2.3. Given e > 0, the subset of described by, 

ll^iXi — X2B1 — Cill < e, 'I 
11 ^ 2 X 3 — X 2 B 2 — C 2 \\ < £, > , 

Xj - X* = Q,l<j <m ] 

where Ai,A 2 ,Bi,B 2 ,Ci and C 3 are some fixed but arbitrary matrix contractions 
in Ain, is a semialgebraic matrix variety. 

Remark 2.1. Matrix sets of the form\2fM and\2f^ provide a connection between 
matrix eguations on words (in the sense of mm) and algebraic/semialgebraic ma¬ 
trix varieties. From this perspective, the computation/refinement of matrix words 
corresponding to the numerical solution of matrix equations on words can he inter¬ 
preted as a discrete analogy of local matrix homotopies like the ones constructed in 
the proofs of the results in 0 

Definition 2.3 (Curved and Flat matrix paths). Given any three hermitian ma¬ 
trices — In < Hi, H 2 , H 3 < In and a function f € C'^([—1,1]), and given any four 
normal contractions Di,..., D 4 in Mn, with D 2 = Ad[e'^*'^^](Zli), D 3 = f{H 2 ) 
and D 4 = f{H 3 ). Let us consider the paths Z{t) := Ad[e’^*‘^'^](Zli) and V{t) := 
f{tH 3 + (1 — t)H 2 ). We will say that Z is a curved interpolating path for Di, D 2 
and we will say that the path V is a flat interpolating path for 0^,0 4 . 


( 2 . 3 ) 


Z" := l{Xi,X2,X3)&Mf, 


Definition 2.4 (@ operation). Given two matrix paths X,Y & C{[Q,V\,Mn) we 
write X ®Y to denote the concatenation of X and Y, which is the matrix path 
defined in terms of X and Y by the expression. 


X®Ys 


X2s, 0 < S < ^, 

I < s < 1. 


Definition 2.5 (Local matrix deformations x v)- Given two matrices x,y € 
Mn we write x y if there is a matrix path z S (^([0,1], M„) such that zq = x and 
= y> */ there is a e-local matrix homotopy X € ^([0, l],Mn) between x and y 
relative to the set S, we will write x -^e,S U will omit the explicit reference to 
S when it is clear from the context. 


It is often convenient to have iV-tuples (or 2A^-tuples) of matrices with real 
spectra. For this purpose we use the following construction, initiated by McIntosh 
and Pryde. If A = (Ai,..., Ajv) is a A-tuple of n by n matrices then we can 
always decompose Xj in the form Xj = Xij + iX 2 j where the Xkj all have real 
spectra. We write 7 r(A) := (An,..., Aiw, A 21 ,..., X 2 n) and call 7 r(A) a partition 
of A. If the Xkj all commute we say that 7 r(A) is a commuting partition, and if the 
Xkj are simultaneously triangularizable 7 r(A) is a triangularizable partition. If the 
Xkj are all semisimple (diagonalizable) then 7 r(A) is called a semisimple partition. 

We say that N normal matrices Ai,..., Ajv G Mn are simultaneously diagonal¬ 
izable if there is a unitary matrix Q € Mn such that Q*XjQ is diagonal for each 
j = \,... ,N. In this case, for 1 < /c < n, let A^^'^{Xj) := {Q*XjQ)kk the {k,k) 
element of Q*XjQ, and set A('')(Ai,..., Ajv) := (A(^^(Ai),... ,A^^'>{Xn)) in C^. 
The set 

A(Ai,...,Aw) := {A«(Ai ,..., Ajv)}i<fc<Ar 
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is called the joint spectrum of Xi,... ,Xn■ We will write A{Xj) to denote the di¬ 
agonal matrix representation of the j-component of A(Xi,..., Xn), in other words 
we will have that 


A.{Xj) = diag 




The following result ([22l L.4.1]) was proved in [22] . 

Lemma 2 . 1 . Given e > 0 there is 6 = -j^s > 0 such that, for any two N-tuples of 
pairwise commuting normal matrices x = {xi, ..., xn) and y := (j/i,..., yjq) such 
that S(x,y) < <5, there is a unitary matrix W such that [Xd\W]{xj),yj\ = 0 and 
max{||Ad[W](xj) — yj\\, ||Ad[lT](xj) — x^jl} < e, for each I < j < N. 


Remark 2.2. The constant Km in the statement of L \2.1\ depends only m. 


3. Local Deformation of Normal Contractions 


Applying functional calculus on commutative C'*-algebras we have that the joint 
spectrum a{X, Y) of any two commuting hermitian matrix contractions X,Y G Mn 
is contained in [—1,1]^, moreover if \\X iY\\ < 1 we will have that a{X,Y) C 
R(0,1) ~ D2. 

Lemma 3.1. Given r G TX and v > 0, there is 6 := S{i',r) > 0 such that for any 
unitary W and any normal contraction D in Mn for n>2, if D — 
diagonal for 2 < r G Z and ai,. .. ,ar G the set {Pj} is an orthogonal partition 
of unity, and aj ^ whenever k ^ j, then there is a unitary matrix Z G Mn such 
that [Z, D] = 0 and ||1„ — WZjf < v whenever ||1L-D1T* — D\\ < 5. 


Proof. Let p := UlTZllT* — D\\ = \\WD — Z1W||. Then there are r continuous 
functions G C(D^) such that Pk := £k{D). By continuity one can find 

S{iy,r) > 0 such that, \\WPkW* - Pk\\ = \\W£kiD)W* - £k{D)\\ = \\£k{WDW*) - 
^k{D)\\ < vl{\f 2 r) < l/(-\/2r), whenever p < S. Since < 1, by perturbation 
theory of C'*-algebras we will have that there is a unitary Wj such that ||1 — 
Wjj| < V^\\WPkW* — Pk\\ < vfr and W*PjWj = WPjW*, and this implies that 
Wj WPj = Pj Wj W for each j. Since {Pj } is an orthogonal partition of unity, we will 
have that W := ^ unitary matrix which satisfies the commutation 

relation \W, D] = 0 together with the normed inequalities. 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


liw-wil 


= WY.P,-Y.W,WP,\\ 

3 3 

= ||^((l„-W,-)WP,)ll 

3 

< ^|| 1 „- 1 T ,|1 

3 

1 A ^ 

3 


— WZ\\ = ||1T — IT|| < V. This completes the 

□ 


Let us set Z := W* then ||1„ 
proof. 
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Definition 3.1 (Pseudospectral retractive approximant). Given (5 > 0 and any 
matrix X G Mn, we say that the matrix X G is a S-Pseudospectral retractive 
approximant fJ-PSRAj of X if ||A — X|| < 5 and a{X) is S-dense in <j{X). 

Lemma 3.2. Given 5 > 0 and any hermitian matrix X G Mn such that ||A|| < 1, 
there is a hermitian S-PSRA Xg of X such that [A, A] = 0„. Moreover, there are 
Ng distinct points {xi ,..., Xn^} o-nd an orthogonal partition of unity {Pi ,..., Pn^} 
such that Xg = 

Proof. Let us suppose that n > |cr(A)| > 2, as the proof is clear for scalar multiples 
of In- Since cr(X) C [—1,1], we can assume for simplicity that Mg := 1 + G 
Z+. Then the finite set Rg{X) := [xk ■= 1)511 < fc < Mg} is 5-dense in 

[—1,1] with |i? 5 (A)| = Mg. Let us set Sg{X) := {ifc := —l+2(fc—2)5|1 < k < Ng-\- 
1} and Pk := X{xk,xk+i]i^) each 1 < fc < Mg, where X(c,d] denotes a continuous 
representation/approximation of the characteristic function X(c,d] of (c, d] on some 
set Sc,d such that {c,d] C Sc,d C R. Then there is Rg{X) := {xj} C Rg{X) 
that is 5-dense in a-{X) with Xj G ixk{j),Xk{j)-i-i] for each Xj G Rg{X) and some 
a;fc(j),ifc(j)+i G S' 5 (A), and there is a set Pg{X) := {Pi,... ,Pns} C M„\{0„} that 
is an orthogonal partition of unity such that [Pj,X] = 0 and for each 1 < fc < 
there is Xj(^k) G Rs^X) such that ||APfe — Xj(^k)Pk\\ < S with Xk(j) x^iy) whenever 
j I. By the previous facts and spectral variation of normal matrices (in the sense 

of [5J11S]) we will have that if we set Xg := J2k Xj{k)Pk with Xji^k) S Rg{X) for each 
fc, then [X,Xg\ = 0„ and dHiX,Xg) < ||A - Xg\\ < maxs, \\XPk - Xj(^k)Pk\\ < S. 
If necessary we can renumber the elements of Rg{X) according to the elements of 
Pg{X). This completes the proof. □ 

Remark 3.1. We will refer to the finite sets Rg{X) and Sg{X) as representation 
and support grids respectively. The finite set Pg{X) will be called a 5-projective 
decomposition for X. 

Lemma 3.3. Given 5 > 0, m > 1 and any m-tuple of pairwise commuting hermit¬ 
ian matrix contractions X = {Xi,... ,Xm) G Mff, there is a m-tuple of pairwise 
commuting hermitian matrix contractions Xg = {Xg^i,... ,Xg^m) € M}f such that 
[Xj,Xg^k] = On and Xgj is a S-PSRA of Xj for each 1 < j,k < m. Moreover, 
there are Ng distinct points {xk,i, ■ • ■ ,Xk,Ns} an orthogonal partition of unity 
{Pi,..., Pns} such that Xg^k = XkjPj for each fc. 

Proof. Since the joint spectrum ct(X) of X is a subset of [—1,1]™ we can apply a 
similar procedure to the one implemented in the proof of L ]3.2l to find for each Xj a 
representation grid Rg{Xj) together with a support grid Sg{Xj) and an associated 
5-projective decomposition Pg{Xj) such that \Pg{Xj)\ = Ngj := \Rg{Xj)\, we 
will have that the set ^^(X) := {Pi,..., PatJ = {Pi,hP 2 ,j 2 ■ ■ ■ Pm,j„, \ Pk,jk £ 
Ps{Xk), 1 < jfc < Ng^k, 1 < fc < rn} is an orthogonal partition of unity, by setting 
(l + 5“^)™ > Ng := |P 5 (X)| and v := l/{Ng — 1), will have that Ng > Ng^k for each 
fc, that V < 5 and that there is a set Rg{Xj) C {xj^k = 2(fc — l)i^ — 1 | 1 < fc < Ng} 
of distinct points that is 5-dense in o-{Xj) for each j and such that for each Xk 
and each Pj G PgfX.), there is Xk,j G Rg{Xj) such that \\XkPj — Xk,jPj\\ < S. 
Moreover, [Xj,Pk] = 0„ for each 1 < j < m and each 1 < k < Ng. By a similar 
argument to the one implemented in the proof of L I3.2I it can be seen that the 


DYNAMICAL DEFORMATION OF TOROIDAL MATRIX VARIETIES 


7 


matrix Xs^k '■= ^k,jPj is a commuting hermitian (5-PSRA of Xj for each j. This 
completes the proof. □ 


Definition 3.2 (Uniform piecewise analytic local connectivity in M™). We will 
say that a matrix variety Z™ C is uniformly locally piecewise analytically 

connected (ULPAC ) if given e > 0, there is <5 > 0 such that for any two m-tuples 
X and Y in Z™ such that S(X, Y) < 6 (X and Y are {6e,d)-close), we have that 
X '^h Y relative to Nq{X, e) n Z™. 

Remark 3.2. It is important to remark that e and 6 must not depend on n. 

Let us define the toroidal matrix varieties whose local connectivity will be studied 
in this document. 


Definition 3.3 (ll"*(n)). The matrix variety D™(n) defined by 



r 


[Xj,Xk]—0n I 

r{n) := < 

(^1,. 

.,Xm)€Mffi 

Xj — Xj = On ,1 < j,k < m '. 

Il^.-ll < 1 J 


will be called the m matrix cube. 

Definition 3.4 (D™(n)). The matrix variety D’"(n) defined by 


D’"(n) := {Z, 


GM" 


[Zj, Zk] — [Zj, Z*] — Or, 


< 1 


,1 < j,k < m 


will be called the m matrix disk. 

Definition 3.5 (T’”(n)). The matrix variety T”^(n) defined by 


T'"(n) 




\Uj, Uk] = o„ 

u,u; = u;u, = i„ 


.1 < j,k < m 


will be called the m matrix torus. 


Remark 3.3. It is important to notice that the components of two m-tuples X = 
(Xi,... ,Xm) and Y = (Yi,... ,Ym) in ll"*(n), D'"(n) or T™(n) need not to satisfy 
the commutation relations [Xj,Yk] = 0„ for each j,k in general. 

Theorem 3.1. The matrix variety ll"*(n) is uniformly locally piecewise analytically 
connected. 


Proof. Without loss of generality we can assume that the components of Y are diag¬ 
onal matrices. As a consequence of L ]2.1l we will have that given i/s > 0, there is <5 := 
> 0 such that for any m-tuple X G D'"(n) that satisfies 6(X, Y) < 6 , there 
exists a unitary matrix W such that [Ad[lU](Xj), Yfi = 0 and max{|| Ad[fU](Xj) — 
YjW, ||Ad[fU](Xj) — Xj\\} < vs, for each I < j < N. Since cr(Xj) C [—1,1] D cr(y,) 
for each j. By joint spectral variation (in the sense of [25]) and by applying L I2.il 
again, we will have that W*h.{Xj)W = Xj and \\W*h.{Xj)W — A(Xj)|| < Km6 
for each j. Then by L I3.3I there is a m-tuple of pairwise commuting hermit¬ 
ian d-PSRA Xi of A(X) := (A(Xi),..., A(Xm)), together with m representa¬ 
tion grids Rs{Xj) := {xkp, ■. ■ ,Xk,Ns} that are d-dense in cr(Xj) = (t(A(Xj)) 
for each j, and a projective decomposition Ps{A{X.)) := {Pi,..., Pns} such that 
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- ^kjPj, 

[Xs,j,A{X 

,)] = On and ||Xj - W*Xs,]W\\ < S. 

±s ■■= 

(W*Xs,iW,. 

.., W*Xs^mW). We will have that. 

(3.5) 

II A. - 

■A. II < 

^Xs,Xs) 

(3.6) 


< 

di±s,X) + g(X, A(X)) + 5(A(X), X 5 

(3.7) 


< 

(i^m + 2)6 


Let US set := |P 5 (Y)|. If = 1, there are four flat hermitian paths := 

X^+tiXs,j-Xj), X^’^it) := Xs,j+tiXs,j-Xsj), X^’^t) ■= Xsj+tiAiXj)-Xs,j) 
and X‘*^'^(t) := A{Xj) + t{Yj — A{Xj) such that the matrix path X^ := [{X^'^ ® 
X'^’^)®X^’^]®X‘^'^ £ C'([0,1], /'"(n)) solves the problem Xj Yj for each j with 
El = {2Km + 3)(5 and this implies that X Y relative to iVg(X,ei) fl I"^(n). If 
iVe > 2 , then we can apply US.ll to find a unitary W := ZW such that W*Xf,^jW = 
W*XsjW = Xs,j and that || 1 „ — IL|| < £2(N^,m,5) < 2 , this implies that one 
can find a hermitian matrix contraction K such that lY = We have that the 

curved path := Ad[e“®^^^“*)^](X 5 j) solves the problem Xs,j Xsj for 

each j with £2 := 2£2{Ns,m,S). Let us set X^ := [{X^'^ ® X'^’^) ® X^'^] ® X'^’^ G 
(^([0,1] with X^'\ and X^'^ defined as before. Then X^ solves the 

problem Xj Yj for each j with £3 := £2 + {K^ + l)d, and this implies that 
X Y relative to IVg(X,£ 3 )nD™(n). Let us set e := max{£i,e 3 }. This completes 
the proof. □ 

Corollary 3.1. The matrix variety D’”(n) is uniformly locally piecewise analyti¬ 
cally connected. 

Proof. Given any two m-tuples Z and S in D'"(n), such that S(Z, S) < r/2 for some 
r > 0, there are two semisimple commuting hermitian partitions 7r(Z) and 7r(S) of 
Z and S respectively, such that S(7r(Z), 7r(S)) < r. By the previously described fact 
the result follows by applying T I3.il to the corresponding semisimple commuting 
hermitian partitions of any two {6^, S)-close m-tuples in D™(n). □ 

Theorem 3.2. The matrix variety T"*(n) is uniformly locally piecewise analytically 
connected. 

Proof. We have that for any two unitary matrices U,V £ Mn such that [U, V] = 0„ 
there are hermitian matrices Tlu, Hy such that the flat path U{f) := 
satisfies the interpolating conditions C/(0) = U, U{1) = V, together with the con¬ 
straints [U{t),U] = [U{t),V] = On and U(t)U(t)* = U{t)*U{f) = 1„ for each 
t £ [0,1]. We also have that for any unitary U and any hermitian contraction K 
the curved path V{t) := Ad[e*'^‘^](17) preserves commutativity and is a unitary for 
each t £ [0,1]. By the previously described facts the result follows by modifying 
the proof of T I3.il conveniently for any two {6^, S)-close m-tuples in T™(n). □ 

4. Hints and Future Directions 

The implications of the main results in ^ in the numerical solution of matrix 
equations on words will be further explored, some applications to the analysis of 
approximation/refinement procedures that can be implemented to solve matrix 
equations on words and eigenvalue problems in the sense of [n n s El m HI] 
will be the subject of further study as well. 
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The approximation and connectivity technology developed in this paper has 
a natural connection to approximate simultaneous diagonalization of m-tuples of 
matrices (in the sense of |24] 1 and to normal matrix approximation of almost normal 
matrices in the sense of [H US]. The development of numerical algorithms to 
perform these tasks will be the subject of future communications. 

The application and extension of the results in lj3] to the study of normal and 
near normal compressions of normal matrices (in the sense of [El [M] and [2 §9 
and §10]) will be the subject of further study as well. 
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